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Abstract 

Magnetic semiconductor quantum dots with a few carriers represent an interesting model system 
where ferromagnetic interactions can be tuned by voltage. By designing the geometry of a doped 
quantum dot, one can tailor the anisotropic quantum states of magnetic polarons. The strong 
anisotropy of magnetic polaron states in disk-like quantum dots with holes comes from the spin 
splitting in the valence band. The binding energy and spontaneous magnetization of quantum 
dots oscillate with the number of particles and reflect the shell structure. Due to the Coulomb 
interaction, the maximum binding energy and spin polarization of magnetic polarons occur in the 
regime of Hund's rule when the total spin of holes in a quantum dot is maximum. With increasing 
number of particles in a quantum dot and for certain orbital configurations, the ferromagnetic state 
becomes especially stable or may have broken symmetry. In quantum dots with a strong ferro- 
magnetic interaction, the ground state can undergo a transition from a magnetic to a nonmagnetic 
state with increasing temperature or decreasing exchange interaction. The characteristic temper- 
ature and fluctuations of magnetic polarons depend on the binding energy and degeneracy of the 
shell. The capacitance spectra of magnetic quantum dots with few particles reveal the formation 
of polaron states. 

PACS numbers: 78.67.Hc, 75.75. +a, 
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I. INTRODUCTION 



Diluted magnetic semiconductors combine high-quality crystal structures with the mag- 
netic properties of impurities 1 and represent an important class of materials for spintronics 
and quantum information 2 . The ferromagnetic ordering in diluted magnetic semiconduc- 
tors can come from the carrier-mediated interaction between magnetic ion a 3 i 4 5 i 6 . Since the 
carrier density in semiconductor field-effect transistors is a voltage-tunable parameter, the 
ferromagnetic state of the impurities coupled to the carriers also becomes controlled by the 
voltage 6 . Voltage-control of the ferromagnetic phase transition has been already demon- 
strated for the Si/Ge and A 3 B 5 material systems 4 . This ability to externally control the 
properties of magnetic crystals with means other than the external magnetic field may have 
important device applications. 

An important feature of modern nanotechnology is the ability to shape semiconductor 
crystals, designing their quantum properties. Quantum confinement of carriers is expected 
to strongly affect the magnetic properties of crystals since the quantum-confined structures 
can be designed in a way to strongly localize carriers near magnetic impurities. The first 
step toward quantum confinement has been made with the quasi-two-dimensional structures 
where unusual ferromagnetic properties have been described&^Sai. For one- dimensional litho- 
graphic structures, it was found that their transport properties are controlled by the domain 
walls 10 . In parallel, self-organization growth technology suggests zero-dimensional nano-size 
quantum dots (QDs) 11 ^ 2 which can locally store carriers. Moreover, it has been demon- 
strated in many experiment a 13 i 14 i 15 i 16 that the numbers of carriers and wave function of a 
QD can be changed by a voltage applied to specially-designed metal contacts; this may 
permit manipulation of ferromagnetic states in QDs by voltage. Therefore, the combination 
of semiconductor QDs with magnetic impurities looks particularly interesting; information 
in such magnetic QDs can be stored not only in the number of carriers but also in the form 
of the Mn magnetization. Currently, magnetic QDs are a hot topi o 1 ^ 1 ^ 1 ^ 2 ?! 21 ! 22 ! 23 ! 24 ! 2 ^ 26 . 
One important property of Mn-doped nanostructures is that a single particle (electron or 
hole) can strongly alter the ground state of the system, leading to formation of a magnetic 
polaron (MP ) 9 ' 25 i 26 , In the case of a QD with a single Mn impurity, single carriers lead to 
the formation of hybrid electron-Mn state a 19 ! 21 ' 22 ' 24 . 

In a semi-magnetic QD, a localized MP is formed due to the exchange interaction between 
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the spins of Mn ions and a carrier trapped in a QD. The MP localized inside a QD resembles a 
localized acceptor-bound exciton in a bulk semiconductor doped by magnetic imp untie d 27 ' 28 . 
However, self-assembled QDs have important differences: (1) a single QD can trap several 
electrons (holes), (2) the number of particles in a QD can be tuned with the voltage applied 
to a metallic contact, and (3) the confining potential of QDs typically is very different from 
the Coulomb potential. 

Here we describe anisotropic MP states of doped QDs of cylindrical symmetry with few 
holes. The binding energy and magnetization of MP states demonstrate oscillations as a 
function of the number of holes due to the shell structure. The maximum binding energy 
and magnetization occur in the regime of the maximum total spin of the hole subsystem 
when Hund's rule is applied. Note that a recent paper— describes only the odd-even parity 
oscillations in magnetic QDs (these oscillations come from the spin susceptibility of a system 
with a discrete spectrum). The enhanced Mn polarization and strong MP binding in a QD 
with cylindrical symmetry described here originate from the symmetry of the system and the 
Coulomb interaction. For QDs with few particles we also predict transitions from magnetic 
to non-magnetic states when the temperature or exchange interaction varies. In addition, 
we focus here on the hole-mediated ferromagnetism in disk-shaped QDs where the MP state 
is strongly anisotropic: the spontaneous magnetization appears preferentially in the growth 
direction. 

Another interesting question related to nano-scale ferromagnetism is how the coupled Mn- 
hole system develops from the MP behavior toward the Zener ferromagnetic phase transition 
regime in the limit of a large number of particles. To address this question, we will also 
consider spin fluctuations in a QD as a function of the hole number and a formal self- 
consistent solution with critical behavior. These results may help to answer the above 
question. The crossover from the MP regime toward the phase transition behavior can also 
be important from the point of view of device applications. For example, it is essential to 
estimate the minimal number of carries needed to achieve a stable ferromagnetic state in a 
single QD. 

The paper is organized as follows: Section II presents a model of a self-assembled QD, 
Sections III and IV describe the anisotropic MP state with one hole, Section V, VI, and VII 
contain the results on few-hole QDs, and Sections VIII, IX, and X discuss critical phenomena, 
fluctuations, and electrical capacitance of magnetic QDs. 
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II. MODEL 



We model the hole-Mn complex in a self-assembled diluted semiconductor QD with the 
following Hamiltonian: 

where R = (r, z) is the radius vector, z is the vertical coordinate, r = (x,y), p is the mo- 
mentum, and j z and S z are the z-components of the hole and Mn momentums, respectively. 
S z = J2i Si )Z 6(R — Hi), where S ijZ and Rj are the spin and position of the z-impurity, respec- 
tively It is convenient to model the in-plane motion of a hole by the parabolic potential and 
the vertical motion with a square well 29 . So, we write: U(R) = u(z) + mu; 2 r 2 /2, where u(z) 
is the z-confinement potential and ujq is the in-plane frequency. The anisotropic exchange 
interaction in eq. ^implies that the QD is disk shaped and the vertical size of QD, L, is much 
smaller than the in-plane wave function dimension, I = i.e. L <C I. Since we 

consider the only heavy- hole states, the exchange interaction (the last term in eq. becomes 
strongly anisotropic 30 . In our model, the light-hole states are assumed to be strongly split 
from the lowest heavy- hole states in the QD. The single-particle spatial hole wave functions 
and their energies in the absence of Mn impurities are given by ij)n,m — fo(z)Xn x ,n y { r ) and 
e„ jm = Tvu}o(n x + % + 1), respectively. Here, the wave function fo(z) corresponds to the 
lowest state in a square potential well in the z-direction, and Xn x ,n y { r ) are the usual wave 
functions of a 2D harmonic oscillator; n x i y \ = 0, 1, 2, ... are the quantum numbers. 

In the spirit of the mean field theory, we can average the operator over the impurity 
positions and write: 



HL = + U(R) - ^x Mn (R)N j z S z (R), (2) 

where % n (r) is the reduced Mn spatial density in the system, Nq is the number of cations 
per unit volume, and ^(R) is the locally averaged Mn spin: 
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FIG. 1: (a) Model of a Mn-doped self-assembled QD. (b) Transistor structure with magnetic QDs; 
the gate voltage controls the number of holes in the QD layer and therefore the magnetic state of 
QDs. Capacitance of this structure can reveal the magnetic state of QDs. 

where S = 5/2, B s is the Brillouin function, j 2 (R) =< x)\jzK R ~ *) > 

is the averaged momentum of the hole at the Mn position, and \^(R,x) > is the wave 
function of a hole which depends on the spatial and spin coordinates, R and \i respectively; 
X = ±3/2; To in eq. El appears due to the anti-ferromagnetic interaction between Mn ions. 

III. MEAN FIELD SOLUTION 

According to eq. El the hole in a magnetic QD moves in the presence of the effective 
spin-dependent potential: 



U eff = U(R) - ^x Mn (K)N j z S z (R). 



(4) 



In a QD with strong spatial confinement, we can neglect the effect of the second term in 
eq. E] on the spatial wave function. At the same time, we should keep it for the spin part of 
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the wave function. The ground state of a magnetic QD with a single hole has a simple form: 



^=^o,o|T>, S z (R) = SB s ( f{ 3 ^ ), ] z (K) = h 0fi (R) 2 , (5) 

where | |> is the hole state with j z = +3/2. Since (3 < for the Mn-hole interaction, 
S z (R) < 0: at low T, the spins of the hole and Mn ions inside the QD are anti-parallel. This 
is a state of magnetic polaron with the energy: Eq(T) = Tiujq — Eb(T). The second term in 
the above equation plays the role of the MP binding energy: 



E b (T) = ^d 3 R&l x Mn (R)N ^SB s (R,T)}. (6) 

In our definition, the binding energy Ef,(T) < 0. The ground state of a MP is two fold 
degenerate since the states j z = ±3/2 have the same binding energy. 
The total Mn polarization is calculated as 



S tot (T) = [ d 3 R[x Mn N SB s ( f /3 ^ )]. (7) 

JR K B {1 + ±0) 

The corrections to the wave function (jSJ) and the energy © can be found by perturbation 
theory, in which SU e ff = — ^XMn(R)N (3/2)S z (R) is a perturbation: 



»(*,x>= £ <Vv " (fl "' ic/ ' //|,/ "-° > ^(fi)iT>. 

(n x ,n y )jt(0,0) 



5E (T) = £ 



< '0O,o(-R)|^e//|'0n ) m > 



|2 



(n x ,n y )^(0,0) e °'° Cn,r 



(8) 



Now we can estimate the precision of the perturbation theory: 5Eo/Eb is about a few % 
for the typical parameters of the problem, which will be specified below. The precision is 
high because of the orthogonality of the spatial harmonic-oscillator functions ( r )- For 
example, the nonzero matrix element < ^2 ) o(i?)|5Z7 e //|V'o,o > a t zero temperature is about 
1 meV and the perturbation-theory parameter < ip n!m (R)\8U e ff\ipo t0 > /(huj ) ~ 0.02. 
With increasing temperature the parameter < ip n ,m(R)\8U e ff \ipo,o > /(^o) decreases and 
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the precision of the perturbation theory becomes further improved. The above estimates 
tell us that the perturbation theory with respect to the potential SU e ff provides us with 
reliable results. We also note that our results are in agreement with ref.— . 

The ferromagnetic state of a QD depends on the spatial distribution of Mn ions, which is 
given by the crystal growth technology. Diffusion of Mn ions during QD growth is driven by 
strain and composition inhomogeneity and can lead to a strongly non-uniform distribution 
of Mn ions in the system. For example, a study of the InMnAs system^ has shown that 
the Mn ions during the growth process mostly substitute for the In atoms inside a QD. 
Here we are going to use a simple model (fig. QJ in which the Mn-doped region forms a disk 
of radius Rmu inside a QD: XMn(p) = % n f° r P < Run and otherwise. Here p is the 
distance to the center of a QD in the 2D plane. It follows from eq. El that the binding energy 
of a MP and the "robustness" of a magnetic state depend on the overlap between the Mn 
distribution XMn(R) an d the wave function of the hole. Since both of them can be tailored 
and controlled by the growth process the ferromagnetic state can be artificially designed. 

Fig. 121 shows the calculated energy of a MP with one hole for different Mn distributions, 
Rmu = 2, 3 nm and oo. The corresponding Mn magnetization is shown in fig. El The 
following parameters of a magnetic semiconductor have been used: XMn — 0.04, (3N = 
—1.3 eV, To = 3.6 K, and Nq = 15 nm" 3 . The QD geometrical parameters were chosen as 
follows: Iq — 4 nm and L = 2.5 nm. The above material parameters represent a CdMnTe 
QD. 

IV. ANISOTROPY OF BINDING ENERGY AND MAGNETIZATION 

The MP state of a hole in a disk-shaped QD is strongly anisotropic due to the valence 
band structure^. This magnetic anisotropy comes from the heavy-light hole splitting in the 
valence band and reveals itself in the last term in the operator (0). The general solution of 
the one-hole problem can be written as 

^ = ^ 0i0 (a|T> +b\i>), (9) 
where |a| 2 + \b\ 2 = 1. Then the MP binding energy takes the form: 
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E b (T) = -f |^ 3 i?[^ (R)a; Mn (R)iV ^(2|a| 2 - ^Bs ^f^f^ \ (10) 

where < \a\ < 1. Figure 0] shows the MP binding energy as a function of \a\. For the cases 
\a\ = 1 and \b\ = 1, the binding energy magnitude is maximum; for \a\ — l/v^2 it equals to 
zero. Therefore, the ground state of a MP corresponds to the pure states | f> or | |>. This 
anisotropy likely plays an important role in optical experiments with excitons trapped in 
semi-magnetic QDs^. In such experiments, an optically-created electron-hole pair rapidly 
relaxes to its ground state resulting in the formation of a MP with the spin parallel (or 
antiparallel) to the growth axis. 

V. FEW PARTICLE STATES 

In the next step, we study the QD magnetization in the presence of few carriers which can 
be loaded from the metal back contact in a voltage-tunable transistor structure (fig.^i). We 
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FIG. 3: Calculated magnetization of Mn ions as a function of temperature for the QD with one 
hole and different Mn distributions; XMn = 0.04 and Rmu = 2,3 nm, and oo. 

will sequentially consider the few first charged states of a QD starting from = 1 (fig. |3J). 
The Coulomb interaction will be treated within perturbation theory which is valid for QDs 
with a strong confinement; namely, we will assume that the quantization energy of QD, 
fkJo, is larger than the characteristic parameter of Coulomb interaction between particles, 
Ecoui = e 2 /(e I). For the QD parameters specified above, we obtain: Hujq ~ 47 meV 
and Ecoui ~ 29 meV. This simplified perturbation approach is very convenient and was 
successfully used for description of experimental data in several publication o 13 ! 15 ! 32 ' 33 . Note 
that the parameters of the QDs studied i n 13 ' 15 are close to those used here. 

In the perturbation approach, we will neglect the Coulomb-induced mixing between shells; 
at the same time, we will calculate exactly the Coulomb-induced mixing within shells by 
diagonalizing the corresponding matrix. Thus, the Coulomb correlations will play a very 
important role for certain states, such as the Hund's states with two particles in the p-shell 
(nh = 4) and three holes in the <i-shell (n=9). 

For some derivations, it will be convenient to treat this problem using the second quanti- 
zation approach. In this approach, the z-component of the total angular momentum of the 
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FIG. 4: Spin anisotropy of the MP energy as a function of the probability to find the hole in the 
state j z = +3/2; XMn = 0.04 and Rmu = °o- 

hole subsystem is given by 



where c+ - z is the creation operator for the single-particle state (7, j z ). Here 7 stands for 
the pair of orbital quantum numbers; 7 may be n x ,n y or n,m z depending on the choice 
of wave functions; here n and m z are the radial quantum number and the orbital angular 
momentum, respectively. For the few-particle wave functions, we will employ the following 
notations: \s h s l ]p lth p ld ;p 2 , h P2,h dij, d lyl ; d 2:h d 2 ,i; d 3jh d 3d ; ... >, where s ta) , p iiTa) , and 
dk,i(i) are the occupation numbers for the s-, p-, and d- states with the corresponding spins. 
These occupation numbers can be either or 1. For the indices i and k, i — 1,2 and 
k = 1,2,3. 

VI. QUANTUM DOT WITH TWO HOLES 

If the QD is occupied by two particles (rih = 2), the spin and spatial variables in the two- 
hole wave function can be separated. Then, for a QD with strong confinement, it becomes 
obvious that the z-component of the total hole spin in the ground state is zero (j z ,tot = 0) 
(fig. and therefore the exchange interaction with the Mn subsystem vanishes. However, 




(11) 



7,J z ==F3/2 
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a) 




FIG. 5: (a) Ground-state configuration of a cylindrically symmetric QD with few holes and without 
Mn ions; the state with nh = 4 is constructed according to Hund's rule. These configurations 
remain ground states if the Mn-hole interaction is weak or the QD confinement is strong enough, 
(b) Ground-state configuration for nt = 9 in the regime of Hund's rule; again this configuration 
represents a ground state if the Mn-hole interaction is not very strong. 

if the confinement is not strong enough, the ground state can change with temperature. 
Figure |Hb shows the hole configuration which can become a ground state if the Mn-hole 
interaction is strong enough. This state has the total spin j z jot — 2 * (3/2) and a non-zero 
binding energy in the presence of the Mn subsystem. The ground state transition occurs 
when 

E 1 -E = \E%^\. (12) 

Here E and Ex are the hole energies of the states (a) and (b) in fig. El The hole energies 
E'o(i) also include the contributions from the Coulomb interaction. E^ h ~ 2 is the binding 
energy for the configuration (b) in fig. 13 

The Hamiltonian of the system is now written as: 
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FIG. 6: Ground-state configurations of a QD with cylindrical symmetry and two holes. The 
configuration (a) is realized in the case of a relatively weak Mn-hole interaction (high temperature 
or small Mn density) and the state (b) corresponds to the system with a strong Mn-hole interaction 
(low temperature or high Mn density). 



H nh =2 =J2& + U ( R ^ ~ ^Mni^NoJ^S^Ri) + U Couh (13) 
i=l,2 '* 

where j Zj tot is the z-component of the total momentum (spin) of holes, Tj is the kinetic energy 
of the i-hole, and Ucoui is the Coulomb interaction. For the Coulomb potential, we use the 
usual formula: e 2 /e\R\ — -R2I, where e is the dielectric constant. 

Regarding the Coulomb interaction in eq. El we will treat it again as a perturbation 
assuming a strongly-confined Q D 13 ' 15 i 32 i 33 . Then, we obtain for the energies of the states: 
E = 2fkJ + Ufj and E 1 = 3huj + U^ r - U^ c , where U^ a , and are the direct and 
exchange Coulomb elements, respectively. The Coulomb matrix elements Jjf^^^ do not 
depend on a particular choice of the single-particle wave functions for the p-shell. For the 
state (b) in fig. |HJ the operator (jllj) includes two terms: 

J z ,tot(R) = ~(I^o,o(#)| 2 ^ s ,t + |Vv(#)| 2 c; t £ p ,t), (14) 
where the p-state orbital wave function ip p can be written as: 

ip p = aip + + pip_; (15) 

here ip + ^(<p) oc e ±(/) and \a\ 2 + |/3| 2 = 1, where r = (p,4>). The indices +(— ) correspond 
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FIG. 7: (a) Orbital anisotropy of the MP energy with = 2 and j z tot = 3. The minimum 
energy corresponds to the state with \a\ 2 = 1/2 in which the angular part of the spatial probability 
distribution is most inhomogeneous, i.e. oc cos{<p + </> ) 2 . (b) The energy of the states 

jz,tot = 3 and j z ,tat = (figs. 6a and b) as a function of temperature; the ground state of the 
system changes with temperature. Rmu = oo and x e ff = 0.04 

to the wave functions with the orbital angular momentum m z = +1(— 1), respectively. The 
binding energy of the MP with two holes 

E b (T,a) = -| Jj 3 R[F 2 (R)x Mn (R)N ^SB s ( ^p^l )], (16) 

where F 2 (R) = |^o,o(-R)| 2 + 1^(^)1 2 i s the particle density in the system. The binding 
energy now depends on the parameter a. We find numerically that the ground state is 
realized for \a\ 2 = 1/2. This ground state corresponds to the wave function with the most 
inhomogeneous p-orbital as a function of <fi: ip p oc cos(<f) + <po), where 0o is an arbitrary 
phase. It can be understood as follows: in the ground state, the spatial wave function of 
holes should be most localized because a strongly localized hole can better control the Mn 
spins and more strongly lower the total energy. This case resembles somewhat a self-trapped 
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MP in the systems with translation invariance^*. We note that the ground state of MP is 
degenerate since the phase 0o is arbitrary in a system with cylindric symmetry. However, if 
the QD confinement is anisotropic, the phase O w iU be fixed by the anisotropy of the QD 
potential. 

In fig. [7k., we show the MP binding energy as a function of |a| 2 . When \a\ 2 = 1 or 0, 
the p-hole state has the orbital angular momentum m z = +1 or m z = — 1, when \a\ 2 = 1/2, 
it has linear polarization. For a QD with the parameters rrihh = O.lmo, fu^o = 47.3 meV, 
and e = 12.5, the Mn-hole system demonstrates a critical temperature at which the ground 
state changes from the state j Ztto t = 3 to the state with j X)tat = (fig. Ujp). This critical 
temperature corresponds to the solution of eq. The existence of this transition depends 
on the particular choice of parameters of a QD. For example, if the density of Mn impurities 
is small enough, the ground state will be always jz,tot = and the transition will not occur 
at all. 

VII. QUANTUM DOT WITH A SMALLER Mn DENSITY AND n h = 1 - 6 

We now consider a QD with rih ranging from 1 to 6. As was pointed out above Coulomb 
correlations will become very important for certain states. In addition, we will assume that 
the Mn density is a few times smaller than in the previous calculations and the ground- 
state transitions of the type shown in fig. do not occur. Therefore, the ground state 
configurations coincide with those in a non-magnetic QD (fig. |3J). The calculated binding 
energies for the lower doping x e ff = 0.01 are shown in fig. [HI 

We now proceed to the case of three holes {rih = 3). In this case, the ground state 
has two holes in the s-shell and one in the p-shell. The closest excited state has one hole 
in the s-shell and two spin-polarized holes in the p-shell (j z jot = 3 * (3/2)). Despite the 
Mn-hole interaction, the configuration j z ^ot = 3 * (3/2) remains an excited state since the 
additional MP binding energy in this state is about 20 meV and significantly lower than 
the quantization energy fvjj for our parameters. In the ground state with = 3, a nonzero 
angular momentum of the hole system comes from the unpaired hole in the p-state. The 
result for the MP binding energy is given by eq. HKlwith Fi(R) = \ip p (R) | 2 . Again the ground 
state is anisotropic due to the degeneracy of the p-shell. 

To describe the ground state of a QD with rih = 4, we should apply Hund's rule. In the 
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FIG. 8: Calculated binding energy of a MP as a function of temperature for various charged states 
of the QD: nh = 1,3,4 and 5. The binding for the states with completely filled shells (n^ = 2,6) 
is zero in our approach; Rmu = re- 
state rih = 4, the s-shell is completely filled, whereas the p-shell is occupied by two holes. 
According to Hund's rule, the ground state of the many-particle system without the Mn 
subsystem should have two particles with parallel spins in the upper shell and the maximum 
spin 3. The corresponding wave function should be found by mixing the Slater determinants 
related to the p-shell and diagonalizing the Coulomb matrix. The states with the smallest 
energy form a triplet: 



II, 1:0,1:1, > +11, 1:1, 0:0,1 > , . 

|1,1;1,0;1,0>, |1,1;0,1;0,1>, 1 ' ' ' ' ' ^ " ( 17 ) 

For the case nu = 4, the first excited states within the p-shell are 1 1, 1; 1, 1; 0, > and 
|1, 1;0, 0; 1, 1 >. The above wave functions should be represented by Slater determinants 
composed of the orbitals tp s and ip±. The state with the largest energy has the configuration: 
|i,i;0,i;i,o>Hi,i;i,0;0,i> j The energies Q f t he above states are given by: E - U^ c , E , and 

E + Upp C , where U^ c is the exchange integral and E involves single-particle energies and 
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some Coulomb interactions. In the absence of the Mn subsystem, the difference of energy 
between the ground state (j z ,tot = 3) and the first excited states (j z ,tot = 0) is equal to 
the exchange energy between p-states: U™ c = {e 2 /e)jdR\dR f +{Rl) ^-f^y^ R2) ^- {Rl) . 
Moreover, the energy of the state = |1, 1; 1,0; 1,0 > will be lowered due to interaction 
with Mn spins and this lowering can be very significant. The reason is the 2- fold increase of 
the total hole spin: 

< *o&,tat(fl)l*o >= 2 l\MR)\ 2 - (18) 

The MP binding energy for this state is 

E b , nh= ,(T) = -f / d 3 R[\MR)\^Mn(B.)3N SB s ( p^^ )}. (19) 

We note that the increase of the binding in eq. EH comes from the factors 2 before and inside 
the Brillouin function. Also, the MP energy (eq. [T9"j) in the case of = 4 does not depend on 
the particular choice of the single particle functions; in the coordinate representation, the p- 
electrons are described with the anti-symmetric wave function: \I/o oc sin(<fii— 02)| ti> | ?2>- 
The magnitude of the calculated MP binding energy demonstrates a strong increase (about 
two times) (fig. EI)- The Mn magnetization also increases for the case nh = 4 (fig. EJ). 
The ground state nn = 5 is again non-uniform as a function of angle 0: 

^ h=5 = a\l, 1; 1, 1; 1, > +/3\1, 1; 1, 0; 1, 1 > . (20) 

When we ignore inter-shell mixing, the results for = 5 become similar to those for nh = 3. 
The magnetization comes from one unpaired p-hole: 

< *o\jz,tot(R)\*o >= +^K+ + /#-| 2 - (21) 

The lowest MP energy is obtained for the "linearly-polarized" state of holes (a = 1/V2) as 
in the case of rtu = 3 (fig. EJ)- 

The state nh = 6 has completely filled s- and p-shells and spontaneous Mn magnetization 
does not appear for QDs with a relatively weak Mn-hole interaction (the inter-shell mixing 
is neglected). 
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FIG. 9: Calculated magnetization of Mn ions as a function of temperature for various charged 
states of QD: = 1,3,4 and 5. The magnetization for the states with completely filled shells 
(n/j = 2,6) is zero in our approach; x e fj = 0.01 and Rmu = °o nm. 

VIII. FORMAL SELF-CONSISTENT SOLUTION FOR A FEW PARTICLE 
QUANTUM DOT: TOWARD THE FERROMAGNETIC PHASE TRANSITION. 

The system of mobile carriers can undergo a ferromagnetic transition if the Mn carrier 
interaction is strong enough and exceeds the anti-ferromagnetic interaction between ions. 
The Curie temperature of Zener ferromagnetism is given by the density of states of mobile 
carriers, the density of Mn ions, and the Mn-hole interaction constant: 

Tcurie oc S(S + l)x ef f N p(E F ) (3 2 / k B , (22) 

where p(Ep) is the density of states at the Fermi level. It is know that the MP state in a 
QD does not undergo the phase transition: the spontaneous polarization in the MP state 
simply decreases with temperature^. This is due to the fact that the system contains just 
one or few carriers. 



17 



We now consider a formal self-consistent solution of the mean field theory in a manner 
similar to the Zener theory. For rih — 1, we employ eq. Eland obtain a self-consistent integral 
equation 



SJR) = SB s ( P/ 3 < Jz ^ ) (23) 
y ' y k B (T + T ) 1 v ; 



where the < j z >t is now averaged over the states of a hole: 



3 6 ^B^ g ^B^ 

< > r = Za(y) , (24) 

where the quantity 

Ai = -f \ J d 3 R[^l x eff N S z (R)} (25) 

plays the role of binding energy for the hole state +3/2; Z s = e k B T + e k s T is the partition 
sum for the s-shell. The inter-shell mixing is again ignored. Assuming a homogeneous spatial 
Mn distribution and integrating eq. 1221 we obtain: 



Ai = -^-Sx eff N Q h{^l (26) 

where /i(Aj) = J B$ ( ^b(t+To) )d 3 R- This simple equation leads to the critical behavior: a 
solution with nonzero magnetization (Ai > 0) exists if T < T crit where 



knT^ 1 = J (| ) 2 jl^P^Xe ff N J ffiRifa - k B T (27) 

This critical behavior occurs in the regime of strong fluctuations in the MP state and 
obviously is incorrect. It is known that MPs do not exhibit a critical behavior. However, 
it seems to be interesting to compute how this formal self-consistent solution develops with 
increasing number of holes in the QD. First, the self-consistent approach becomes more 
reliable with increasing the number of particles. Second, eq. 123 gives a useful estimate of 
the characteristic temperature of spontaneous magnetization for the MP states in a QD. 

For the cases of nh = 4 and nh = 9, we have found similar self-consistent solutions for the 
energy and magnetization by solving the corresponding nonlinear equations. For example, 
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in the case of nh = 9, the ground-state configuration is constructed according to Hund's rule 
(fig. Eb) and the self-consistent problem is reduced to a system of two nonlinear equations 
of two variables: X\ = f d 3 R[ip^ S Z (R)} and x 2 = / d 3 R[ip^ ,S Z (H)], where ^,0 an d V'd.H- 
are the d-orbitals with m z = and +2, respectively. Again we ignored inter-shell mixing 
and diagonalized the Coulomb matrix inside the p- and d-shells. Then, using the obtained 
energies, we constructed the partition sums for the p- and d-shell. In figure ITUb . we show 
by dashed lines the averaged energies of MPs calculated within the self-consistent approach. 
In the same figure, the solutions without the self-consistent averaging of the spin inside the 
Brillouin function are shown as solid curves. 

Equation (|27|) gives a characteristic temperature at which spontaneous Mn polarization 
exists in a QD. From fig. El we see that the critical temperature increases with number 
of holes. However, this increase is not very strong. The reason is that the p- and d-shells 
have several states and some of these states have opposite spins or no spin. These states 
contribute to the partition sum and therefore effectively reduce the magnetization. 

It is interesting to compare the typical temperature of MPs with that of ferromagnetic 
phase transitions in bulk. For CdMnTe, it was suggested that the Curie temperature for the 
highest hole densities can be as high as a few Kelvins^. In a QD with about 20 Mn ions 
considered here (x e // = 0.01, / = 4 nm, and L z = 2.5 nm), the typical temperature given by 
eq. Ellis about 22 K, an order of magnitude larger than that in bulk. This is also consistent 
with previous papers on MP o 19l2 ° . We can also note the important differences between the 
above equations for Tcurie and T crit . The critical temperature is proportional to the first 
power of the interaction (3 and depends on the QD localization length (T crit ~ Z -3 / 2 ). 



IX. FLUCTUATIONS AS A FUNCTION OF n h 

To better understand the behavior of MPs at high temperature, we now compute the 
relative fluctuations of the hole subsystem: 



v < (jz,tot- < jz,tot >T) 2 >T t 

Sj z ,tot(T, n h ) = -5! , , (28) 

< 3z,tot >T 

where < ... >t means thermal averaging over the many-particle states assuming a given Mn 
spin distribution. For the Mn spin distribution we will use eq. El calculated for a certain 
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FIG. 10: (a) Calculated average spin of the hole subsystem (solid curves) and its relative fluctu- 
ations (dashed lines) for the states rih = 1,4, and 9. One can see that the ferromagnetic state 
becomes more stable with increasing number of holes, (b) The solid curves show the calculated 
binding energy of a MP as a function of temperature for three charged states of a QD (n>h = 1,3, 
and 9) with maximum binding energy; XMn = 0.01 and Rmu = oo. The dashed curves show the 
results of self-consistent mean-field theory which demonstrates critical temperature. The dashed 
straight line is the thermal energy. 

quantum state of the hole subsystem. 

In the low-temperature regime, the relative fluctuations are much weaker than the average 
spin of holes (fig. ITTij) . With increasing temperature, the fluctuations grow and exceed 
< jz,tot >t- Figure [TUl shows clearly the tendency of stabilization of the MP state with 
increasing number of holes for the states = 1, 3, 9. This stabilization is consistent with 
the increase of binding. Clearly, methods beyond the mean-field theory are required to 
understand the mechanism of crossover from the MP behavior to the Zener ferromagnetic 
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FIG. 11: (a) Calculated energies of different charged states in the layer of semi- magnetic quantum 
dots with XMn = 0.01; Rmu = °o. The labels show the charged states of the system. 

phase transitio n 28 ! 36 . 

X. CAPACITANCE OF QD SYSTEMS WITH MAGNETIC POLARONS 

One efficient method to study quantum states of QDs is capacitance spectroscopy. Such 
spectroscopy is typically performed at a nonzero frequency uo. The capacitance of the QD 
structure (fig. ^d) includes the contribution of charges trapped inside the QD layer—: 



where V g is the gate voltage and D(V g ) is the effective density of states in the QD layer. 
The latter is denned as 



6C = e 2 D(V g ), 



(29) 



D(V g ) = N dot 



dn h 



(30) 
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FIG. 12: (a) Calculated capacitance spectra for non-magnetic (a) and magnetic (b) QDs; x e fj = 
0.01. The broadening of the peaks was taken as 16 meV. 

where n h is the number of holes trapped in a single QD and N dot is the 2D density of QDs. 
This approach to the capacitance is valid at low frequencies lot <C 1, where r represents 
both the tunneling time and the relaxation time to form the ground MP state. In other 
words, this approach assumes that the weakly coupled system " QD + metal contact" has a 
short relaxation time and always remains in its ground state while the gate voltage changes 
in time as V g + 5V g cos{ojt) (V g ^> SV g ). 

According to the simple model of a field-effect structure with a QD layer—, the energy 
of a single particle in the QD is written as: E (V g ) = E sp + |e|7V^, where E sp is the single- 
particle energy inside the QD, 7 = d\jdi is the lever arm coefficient, and d\ and g?2 are the 
dimensions of the structure (fig. Q). For example, 7 ~ 1/6 in ref.— . For convenience, we 
assume that the Fermi energy of the metal contact is zero and loading of the first hole to 
the QD occurs at zero bias. Then, the first charged states of the QD (n^ = 1 and 2) have 
the energies: 
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E 1 (V 9 ) = 1 \e\V g + E^-\ 
E 2 (y a ) = 2i\e\V a + U?. (31) 



The energies of the next charged states (n^ = 3 and 4) are 



E 3 (V g ) = hw + 3 7 |e|V^ + U coul , 3 + E^~\ 

E 4 (V g ) = 2frw + 4 7 |e|V fl + U coulA + E™ h= \ (32) 

The Coulomb energies Uc ou i,n h m the above equations should be calculated for the ground- 
state configurations shown in fig. EK- Figure [TT] shows the calculated energies of the first 
changed states related to the s- and p-shells for non-magnetic QDs. As was realized in several 
experiments, the ground state of the system changes with voltage: a QD sequentially traps 
1, 2, 3, ... particle a 13 i 15 . Then, the quantity D(V g ) and the capacitance demonstrate peaks at 
the voltages of the ground-state transitions. In real QD systems these peaks are broadened 
due to the nonzero size dispersion in a QD ensemble. The effects of the ferromagnetic 
interaction are clearly seen in the calculated capacitance spectra (fig. H2J). The spacing 
between the p-orbital peaks rih = 4 and 5 becomes strongly increased. It comes from the 
strong ferromagnetic coupling in the regime of Hund's rule for the state rih = 4. The spacing 
between the two s-orbital peaks is also increased due to the MP effect. At the same time, 
the voltage interval between the s- and p-related structures becomes reduced. This is again 
due to the exchange interaction. 

According to fig. the magnetic and non-magnetic QDs with the parameters chosen in 
this paper show the same order of peaks in the capacitance spectra. With increasing Mn- 
hole interaction, the situation can change and the sequence of peaks can became different 
for the magnetic and non- magnetic systems. The reason is that, in magnetic QDs, more 
carriers can be trapped to achieve the minimum energy. 

In addition to the characteristic behavior of the inter-peak spacings, the MP effect in the 
capacitance spectra can be recognized by varying temperature or by applying an external 
magnetic field. With temperature, the peaks related to the most bound MP states (rih = 
1,3,4,5) will have the strongest temperature dispersion. As for the magnetic fields, one 
possibility is to use the in-plane field and to suppress the Mn-hole interaction. Then, the 
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characteristic magnetic dispersion of the peak positions will reveal the MP binding energy. 
The behavior of the peaks in a perpendicular magnetic field can also be revealing as was 
shown in other publications^. 

XI. DISCUSSION 

The model with exchange interaction of the type j z S z ^5(Rh — Rl) assumes that the 
magnetic impurity does not create any spin-independent potential. This model is widely 
applied for II- VI semiconductors where individual impurities do not form bound acceptor 
states. This is in contrast to the GaAs system. In the GaAs crystal, a magnetic Mn impurity 
forms a deep acceptor state (about 110 meV above the top of the valence band) and therefore 
the model of impurity in bulk GaAs should incorporate the effect of the spin-independent 
attracting potential. Treating the typical /// — V QDs realized in the InGaAs system we 
should take into account two factors: (1) the possible lower binding energy in the InGaAs 
system for the Mn acceptor stated and (2) the spatial confinement in a QD. To understand 
the importance of the acceptor potential, we should compare the localization length of the 
Mn-acceptor in bulk with the size of the QD. If the QD dimension is smaller than the 
Mn-acceptor size, the QD can be treated without the acceptor potential. Simultaneously, 
the spin-dependent exchange interaction should remain in the model since it leads to the 
formation of the MP state. In the opposite limit of a weak QD confinement, the QD potential 
can be treated as a perturbation; this case was recently analyzed in ref.— . To summarize, 
if the electronic size of a QD becomes smaller than the dimension of the acceptor states of 
the Mn impurities inside a QD, the simple model of contact exchange interaction becomes 
applicable. This suggests that, under certain conditions, InGaAs QDs can also be treated 
with the simple model used in this paper. 

To conclude, we have calculated the MP energies and associated capacitance spectra of 
QDs with a few holes in the presence of the Mn-hole exchange interaction. The system 
studied in this paper exhibits several features coming from the joint action of the Mn-hole 
exchange coupling, quantum confinement, and Coulomb interaction. 
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